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We study a model thin film containing diluted bilayer structure with the RKKY long-range 
interaction. The magnetic subsystem is composed of two magnetically doped layers, separated by 
an undoped nonmagnetic spacer and placed inside a wider film modelled by a quantum well of 
infinite depth. We focus our study on the range of parameters for which the antiferromagnetic 
coupling between the magnetic layers can be expected. The critical temperatures for such system 
are found and their dependence on magnetic layer thickness and charge carriers concentration is 
discussed. The magnetization distribution within each magnetic layer is calculated as a function of 
layer thickness. The external field required to switch the mutual orientation of layer magnetizations 
from antiferromagnetic to ferromagnetic state is also discussed. 
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C/5 



C 

O 

o 



> 

oo 
oo 

o 
o 

On 
O 



X 



I. INTRODUCTION 

The phenomenon of mutual coupling of magnetic lay- 
ers separated by a nonmagnetic medium has been of 
great interest since its discovery in Fe/Cr structure^ 
both due to its immense application potential and the 
highly nontrivial physical backgroundi^^^ 7 ^. In par- 
allel to studies of metallic systems, the progress in the 
field of thin films made of diluted magnetic semiconduc- 
tors (DMS), such as the most representative (Ga,Mn)As, 
offers some new possibilities to study the diluted mag- 
netic systems with confined geometry. Such materials 
are promising from the future spintronics point of view^ 
since they can assure integrability of novel spin-current 
devices with conventional structures of semiconductor- 
based electronics J£ In DMS an indirect coupling be- 
tween localized spins mediated by charge carriers is of 
paramount importance for the possibility of magnetic or- 
dering. The Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction mechanism was suggested to describe the 
properties of quantum- well based DMS systems a decade 
agoAi 

The problem of interlayer coupling in multilayers and 
superlattices of DMS has been subject to various theo- 
retical studies. Giddings et alr^ performed some calcu- 
lations within a single parabolic band approximation for 
kp method and local density approximation to study the 
possibility of the appearence of antiferromagnetic (AF) 
interlayer coupling in the superlattices of (Ga,Mn)As. 
They concluded that the minimization of spacer thick- 
ness accompanied with relatively short superlattice pe- 
riod is advantageous for a strong coupling. This prob- 
lem has also been investigated by Sankowski and Kac- 
man,— who found in their tight binding-based study a 
relative insensitivity of the coupling energy to the mag- 
netic layer thickness. A noticeable range of parameters 
favourizing the AF interlayer interaction has been identi- 
fied. The self-consistent mean-field calculations for mul- 
tilayers were also performed by Jungwirth et al^ 



Numerous experimental works evidenced the fact that 
two magnetically doped layers can be indirectly coupled 
for some interlayer distances. For example, such a cou- 
pling has been observed in (Ga,Mn)As/GaAs trilayers by 
means of magnetotransport and magnetization process 
studies by Akiba et al^, where the quantitative analysis 
was performed yielding the ferromagnetic (F) interlayer 
couplings close to 0.1 mJ/m 2 for the nonmagnetic spacer 
up to 10 ML thick. The superlattice structures of the 
same DMS have been investigated by Mathieu et aL— , 
where the critical temperature was found to oscillate with 
the spacer thickness (up to 9 ML) and thus the presence 
of interlayer coupling was deduced. Short-period super- 
lattices formed the basis of neutron reflectivity studies of 
Szuszkiewicz et aliU, who confirmed the ferromagnetic- 
only interlayer coupling (mainly for 6 ML spacer) . In the 
work of Kirby et aL— neutron reflectometry and magneti- 
zation curve studies of trilayer structures were performed, 
in which the magnetization of top and bottom layer has 
been resolved. This indicated distinct ferromagnetic cou- 
pling through a spacer of large thickness (about 20 ML). 
It was also excluded that the spacer has been plagued by 
diffusion of the magnetic impurities leading to direct cou- 
pling mechanism. On the other hand, the results of Ge et 
al^ show the difference of critical temperatures for top 
and bottom layer in a trilayer system, a difference which 
had a tendency to vanish when the nonmagnetic spacer 
thickness was reduced (and a single common transition 
temperature has been observed for the spacer of about 
10 ML thick). Let us note that in that paper it has been 
stated that the sign of coupling energy might be spa- 
tially inhomogenous - either F or AF for various areas on 
the surface of the sample, which was deduced from the 
observation of planar Hall effect. However, it was not 
until last year that the AF coupling between layers in 
(Ga,Mn)As/GaAs:Be was clearly observed by means of 
DC magnetometry as well as neutron reflectivity studies 
performed by Chung et al^ The critical temperature of 
the system was about 50 K. All the discussed interlayer 
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FIG. 1: Schematic cross-section of an ultrathin film contain- 
ing two magnetically doped layers. In this figure the total 
thickness of the film amounts n=ll, whereas the thickness 
of each magnetic layer is nz,=3 and the thickness of non- 
magnetic spacer is ns=3. 



coupling signatures have been assigned to the carrier- 
mediated interaction mechanism. 

The above results encourage for further search for the 
possibility of obtaining AF interlayer carrier-mediated 
coupling. Therefore, it seems purposeful to devote a 
theoretical study to the trilayer structure, which can be 
formed by two magnetic layers immersed in the thin film 
and separated by a nonmagnetic spacer, with the aim 
to discover some unique properties of RKKY interaction 
in the ultrathin film. From the practical point of view 
it is interesting to find a range of carrier concentrations 
and geometric characteristics with AF coupling. In par- 
ticular, the influence of the magnetic layer thickness on 
the coupling energy and the critical temperature of the 
system poses an important question. In order to have a 
deeper insight, also the spatial distribution of magneti- 
zation within the magnetic layers appears worth-while to 
examine. 

For the above reasons, in this paper our purpose is to 
study a model trilayer system with the RKKY interac- 
tion. The theoretical model is based on the quantum- 
well approach to the thin film for which the RKKY ex- 
change coupling has been re-derived. The numerical cal- 
culations have been performed for some realistic material 
constants, typical of DMS systems. 



II. THEORETICAL MODEL 

We consider a model thin film consisting of n mono- 
layers, where some of the planes can be doped with mag- 
netic impurity ions. The doped planes form two magnetic 
layers, separated by a non-magnetic spacer (Fig.l). Each 
magnetic plane consists of ul doped atomic planes, while 
the undoped spacer has the thickness of ns monolayers. 
As a consequence, the inequality 2n^ + ns < n must 



hold. The film possesses fee structure and its surface 
orientation is (001). The localized spins situated in the 
lattice site (r.;, z v ) can be shortly denoted by S^, where 
v is one of the atomic planes whose position is given by 
z v = [y— l/2)d, and is the vector in the plane, d is the 
thickness of a single monolayer, and the total thickness 
of the film is D = nd. The choice of structure is moti- 
vated by the fact that the ultrathin films of (Ga,Mn)As 
are most often grown with this surface orientation, on 
(OOl)GaAs substrate. Our film constitutes a quantum 
well of infinite depth in z direction for the charge carri- 
ers. 

For the free-electron model, the perturbation Hamil- 
tonian describing the exchange interaction between the 
single localized magnetic moment at (rj,z„) and the free 
charge carrier at (r, z) is assumed to be in the following 
form: 



H 



( r ,z), 



-Ap{v,z) S* u . 
where A is the exchange constant, 

1 -((r-n) 2 +(^-^) 2 )/2 



p(r,z) = 



a 3 (2tt) 
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(1) 



(2) 



and s z (r, z) is the z-component of the free-electron spin 
in the (r, z) point of the film. 

We emphasise the fact that the selected " diffused" con- 
tact potential p (r, z) differs from the usual Dirac delta- 
function, since it has a non-zero dispersion a 1 . The con- 
tact potential frequently used in derivation of RKKY in- 
direct exchange integral is reproduced in the limit a —> 0. 
Such a form of the formula (J5J) is further substantiated by 
the fact that the localized magnetic moments originate 
from the ci-type electronic orbitals which possess some 
spatial extension. Note that such a potential has already 
been applied in studies of bulk DMS.— i 2 ^ 3 . 

The influence of "diffusion" of the contact potential 
on the resulting RKKY coupling integral in ID, 2D and 
3D systems has recently been studied in Ref. [24|. The 
studies resulted in obtaining a finite value of RKKY in- 
tegral for r = in 2D and 3D case. The importance of 
this result has first been noticed for discrete 3D systems 
since the presence of divergency at r — > leads to some 
divergent self-energy corrections , 2 ™ However, it must be 
emphasized that in studies of ultrathin films the accep- 
tance of a finite a value is necessary, even if we do not 
require the self-energy to be finite. The reason is that 
the RKKY exchange integral derived on the basis of the 
contact potential with a = is divergent for arbitrary 
z when r = 0, i.e., in the direction perpendicular to the 
film plane, even though the distance between the inter- 
acting spins is nonzero. This fact precludes the further 
studies of the total energy of the discrete system if a = 0. 

By applying the second-order perturbation calculus we 
obtain the RKKY exchange integral between the local- 
ized spins separated by the distance Ar = — r j | in 
the plane of the film, and Az — \z v — z^\ in the direction 
perpendicular to the film surface. The integral has the 
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form of: 
J (Ar, Az) 



= C 
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In the above formula, <p (I) is the Fourier transform of the 
z-dependent part in the interaction potential p (r, z), and 
it can be presented as <f> (0 = ex P (i^la) for a/d < 
The symbol xi,y denotes the paramagnetic susceptibility 
of the electrons in the ultrathin film, which has been 
derived elsewhere ; 26 ! 27 while xo is the Pauli susceptibility 
of 2D system. The energy constant C in Eq.(3) is given 
by C = 2mA 2 /irh 2 a 4 , where a is the lattice constant. 

The Hamiltonian for the system under consideration is 
assumed to be of Ising-type with the long-range exchange 
interaction J (Ar, Az): 



K = - £ J(Ar,Az)^a,^5^ 



(4) 
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The above Hamiltonian includes the Zeeman term for the 
system embedded in an uniform external field B. The 
localized spins reveal the effective gyromagnetic factor 
9cS — 9s + ge(jtAmT-p /h?D), 2 ~ where gs and g e are the 
gyromagnetic factors of localized spins without RKKY 
interaction and the free-electrons, respectively, rp is the 
number of subbands, splitted due to the discretization of 
Fermi surface in thin film^ 

In order to describe a site-dilution in the system we 
make use of the Edwards operators which possess 

two eigenvalues: 1 when the site is occupied by a 

magnetic impurity, or otherwise. In this paper, for the 
simplicity, we will neglect the correlations of opera- 
tors, i.e., the Warren-Cowley short-range-order parame- 
ter— is assumed here to be 0. It means that the occupa- 
tion of lattice sites by the impurity ions in doped planes 
is completely random. As a consequence, the configu- 
rational averages read (&,^) r = x and (6, M £j>) r = x 2 , 
if the sites (i,fi) and (j, /x) belong to the magnetically 
doped atomic planes and vanish otherwise, whereas x is 
the magnetic dopant concentration. In further consider- 
ations it will be assumed that the indices \i and v run 
only over the magnetically doped atomic planes. 

By taking thermodynamical and configurational aver- 
age of the Hamiltonian ^ for the temperature T — > 0, 
we obtain the ground-state enthalpy H =< Ti > of the 
magnetic system under consideration. The enthalpy can 
be conveniently written as a sum of contributions origi- 
nating from the energy of spin interactions within a sin- 
gle magnetic layer, the energy of interactions between 
two layers and the interaction with an external field. 
For the temperature T — > the magnetization in each 
atomic plane v constituting a specified magnetic layer is 



the same, for instance, \S?) = 171% = ±S for the first 
layer and (Sf „) = vn-i = ±5 for the second one, i.e., both 
layers are fully polarized (ferro- or antiferromagnetically 
each to the other). As a consequence, the total enthalpy 
per one site in a doped atomic plane reads: 



H E\ 



intra . -^intcr 



N N N 



xn L g e sl^BB (mi + m 2 ) . (5) 



Since each layer is ferromagnetically ordered inside, the 
energies of intralayer and interlayer interactions can be 
written in the form of: 



N 



kj centra 



and 
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X mim 2 Winter- 



(6) 



(7) 



In the above equations we denoted the appropriate lattice 
sums by: 



n L — 1 



£intra — 



and 



Winter = n L £ {n L + n S ) 



nr — 1 



^2 (n L - v) [£ (n L +n s + v) + £ (n L + n s - f)],(9) 



u=l 

where 



£^) = ^l iui) j(f^\.d). 



(10) 



The logic index V (w) stands for the parity of lu, which 
refers to the fact that in (OOl)-oriented fec-based thin 
film, the subsequent atomic planes are mutually shifted 
parallelly to the surface. This requires two sets of co- 
ordination zones radii rk, and co-ordination numbers 
Zk, to be used: one if both spins are situated in the 
same plane (or in the planes seperated by an even num- 
ber of interplanar distances), V(u>) — "even", and the 
other for the situation when both spins lie in the planes 
seperated by an odd number of interplanar distances, 
"P(o;) = "odd". Both sets of numbers have to be de- 
termined numerically for the assumed crystallographic 
structure. V (0) requires the point (r = 0, z = 0) to be 
excluded from summation. 

In the system considered, the area per one site in the 
(001) atomic plane is S = a 2 /2 so that the interlayer cou- 
pling energy per unit area of the film is 2x 2 S 2 £i n ter / a 2 ■ 

Let us note that the interlayer coupling energy is often 
defined as the difference in energies of the system for fer- 
romagnetic (E^) and antiferromagnetic (E^) alignment 
of magnetic moments in these layers. In our notation 
such a difference is (£ TT - £ u ) /N = ~2x 2 S 2 £ inter . 
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When we consider the magnetic monoatomic layers 
with the thickness = 1, the expressions © and © 
reduce only to the first terms, without summation over 
v£k On this basis, as a special case, we can study the 
interplanar coupling energy for two magnetically doped 
planes as a function of their separation for different total 
thicknesses of the film. 

In general, for non-zero temperature, the magnetiza- 
tions per one magnetic impurity in each magnetically 
doped plane are not equal for all /i, so that some spatial 
magnetization profile is present in the magnetic layers. 
In the molecular field approximation (MFA), these mag- 
netizations can be obtained as the solutions to the set of 
2til coupled self-consistent equations 



SB S 



(11) 



It should be kept in mind that the index /i here runs 
only over the magnetically doped monolayers. Bs (x) is 
the Brillouin function for spin S and the molecular field 
acting on a given spin in the layer \x is 



x^2,£ - v\)m l 



(12) 



The summation over v in Eq.(12) is restricted to the 
magnetically doped atomic planes. 

From linearization of the above set of equations for 
magnetization, in vicinity of the critical point, we obtain 
the critical temperature of a second-order phase transi- 
tion. It is given by the largest real root of the equation: 



det J = 0, 



(13) 



where the matrix J is of size 2ul X 2nr, , and its elements 
are 



s{s + r 

3k B T c 



£{\n-v\). 



(14) 



Due to the long-range summation occurring in 
£ (|/i — v\) (as seen from Eq.(10)), the critical tempera- 
ture can be calculated from the equation (fl3|) only numer- 
ically. The formula given above can be approximated if 
we assume an uniform magnetization distributions within 
each of the magnetic layers, i.e., when we substitute mi 
for the average magnetization in each atomic plane inside 
the magnetic layer 1 and mi for the layer 2, respectively. 
Then the critical temperature can be calculated from the 
approximate formula 



S(S+1) 2g intra - n L £ (0) ± £ u 
3 til 



(15) 



where the "+" sign corresponds to the Curie temper- 
ature and is valid for ferromagnetic interlayer coupling 
(Winter < 0), while the "— " sign is for the Neel temper- 
ature and is valid when the magnetic layers are coupled 
antiferromagnetically (Sinter > 0). The highest accuracy 



of the uniform approximation (15) is for small magnetic 
layer thicknesses n^. 

When AF coupling exists between two magnetic lay- 
ers, the mutual orientation of their magnetizations can be 
switched from antiparallel (antiferromagnetic) to paral- 
lel (ferromagnetic) one by applying an external magnetic 
field along the direction of their magnetizations. This 
kind of field-induced phase transition can take place be- 
low the Neel temperature. The critical field H c required 
to force the reversal of magnetization at T — » can be 
determined by equating the total enthalpies of the sys- 
tem with FM and AFM orientation of magnetizations in 
the presence of that external field. The enthalpies are ob- 
tained from the formula ([5]) , which leads to the condition 
for H r : 



Mo 



g e gfj, B SNxn L 



(16) 



The above expression can be written in the form more 
convenient for numerical calculations: 



where 



Ho He — Bo 



B = 



x \£ u 



n L C 
CS 



(17) 



(18) 



gosfJ-B 

is a material-dependent constant. 

III. NUMERICAL RESULTS AND DISCUSSION 

Since the existence and characteristic of RKKY inter- 
action essentially depend on the electronic structure of 
the system under consideration, we make it a point of 
a brief discussion for the ultrathin film. Confining the 
charge carriers in a quantum well results in discretiza- 
tion of the Fermi surface, which for such a system is 
composed of a finite number of circular 'slices', each of 
them corresponding to one 2D energy subband. The de- 
tailed description of the electronic structure of an ultra- 
thin film is presented in Ref. [2^. In the Fig. [5J the crucial 
quantity for studies of electronic properties, namely the 
density of states (DOS) at Fermi level is plotted against 
the number of monolayers which make up the whole film. 
The method of calculating DOS for an ultrathin film has 
been presented in details in Ref. 123 . In further studies we 
make use of a normalized charge carriers concentration 
p = n r d 3 . In the main plot, it is assumed that the number 
of charge carriers in the film is fixed and their normalized 
concentration equals p = x/n, thus carrier concentration 
decreases when the film is made thicker. This assump- 
tion corresponds to the situation in which only two out of 
n atomic planes contain impurities which serve as charge 
donors, and the concentration of impurities within each 
doped plane is equal to x. Three impurity concentrations 
were selected, namely x =0.025, 0.050 and 0.075 all be- 
longing to the physically relevant range in DMS. For the 
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FIG. 2: Density of states at the Fermi level in an ultrathin 
film, normalized to the bulk value, as a function of number 
of monolayers . The charge carriers concentration is assumed 
p = x/n for three representative values of x. In the inset, 
a similar plot is shown for the charge carrier concentration 
p = x/2. 



inset we chose constant p — x/2, which is valid when all 
the atomic planes are doped with concentration x. For 
convenience, the values of DOS were normalized to the 
corresponding values for the bulk case (for the appropri- 
ate value of p). 

It is visible that DOS undergoes discontinuous jumps at 
some values of n. This is clear manifestation of quan- 
tum size effects (QSE) as each of the jumps occurs for 
the charge carrier concentration at which the next energy 
subband becomes occupied by the carriers, starting from 
a single 2D-like subband for the lowest concentrations. 
Although p is inversely proportional to the film thickness, 
the number of occupied energy subbands increases while 
increasing n. The distance between subsequent jumps be- 
comes shorter by increasing the concentration x of charge 
carrier donors. Between the jumps, normalized DOS de- 
creases as n~ 2 / 3 , which reflects the fact that for a thin 
film the total DOS is a sum of (equal) contributions from 
each occupied 2D-like subband. This 2D-like DOS per 
one charge carrier is energy- independent, so that dividing 
by the 3D DOS per one charge carrier (proportional to 



P 



-2/3 



) yields this behaviour. The dependence presented 
in the main plot in Fig. [2] is a result of interplay between 
the changes caused solely by the change in quantum well 
width and by decrease of the carrier concentration when 
the film becomes thicker. In the inset, for constant car- 
rier concentration, we observe that the magnitude of QSE 
tends to vanish more rapidly for thicker films. For the 
film thicknesses n for which QSE manifest, we can expect 
similar discontinuous bahaviour of other physical prop- 
erties which depend crucially on the DOS at the Fermi 
level. 

Having a special case of the formula © for til = 1 (as 
in Ref 251), we can study the interplanar coupling ener- 
gies £ (Az/d). In the formula (jTOJ) , the summation of ex- 
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FIG. 3: Interplanar coupling energy for an ultrathin film con- 
taining n monolayers, as a function of separation between the 
planes, for the charge carriers concentration p = x/n and 
x = 0.05. 



change integrals over the in-plane coordination zones has 
to be performed. Some cut-off radius must be assumed, 
large enough to assure satisfactory convergency of £. It 
should also be emphasized that calculating each value of 
exchange integral J (rj, , tod) constitutes a remarkable nu- 
merical task itself. Therefore, a noticeable computational 
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FIG. 4: Interplanar coupling energy in an ultrathin film for 
7i = 20, x = 0.05, and the charge carriers concentration 
p — x/n. Full points denote the values obtained from sum- 
mation of the RKKY thin film-modified exchange integral. 
The empty points represent the values obtained from sum- 
mation of the usual RKKY bulk formula. In the inset, the 
very exchange integrals J (r = 0, z) are plotted - the solid line 
is for the RKKY interaction in thin film, while the dashed 
line corresponds to the bulk 3D RKKY coupling. 
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FIG. 5: Interplanar coupling energy for an ultrathin film char- 
acterised by n — 20, x — 0.05, and for three representative 
charge carrier concentrations. 



effort is required to study the inter- and intralayer inter- 
action energies so that we performed our calculations on 
the multi-CPU cluster. We performed the summation 
up to k =1000 coordination zones (which corresponds to 
the cut-off radius r/d=82.0 for even interplanar distances 
and r/d=86.8 for odd interplanar distances). In addition, 
we applied some averaging intended to remove the com- 
ponent oscillatory in cut-off radius. Such a procedure 
provides sufficient convergency of the resulting energies. 
Let us emphasize that accepting a finite spatial exten- 
sion of an interaction potential ([T]) is necessary to obtain 
finite interplanar coupling since the formulas contain the 



exchange integral J (0, z) which would be divergent oth- 
erwise. In further calculations we assumed cx/(i=0.35, 
which could correspond to 1 A for (Ga,Mn)As (which 
we assume to be a realistic estimate of the spatial exten- 
sion of d orbital carrying the localized spin). 

It is interesting to follow the evolution of the energy 
dependence on interplanar distance when the film thick- 
ness n increases. We make an assumption that the charge 
carrier concentration is p = x/n. In Fig. [3] we present 
£ (A/d) as dependent on interplanar distance for vari- 
ous n and for x — 0.05. It is remarkable that all the 
curves are symmetric with respect to the centre of the 
film, i.e. £ (u>) = £ (n — lu). It is a result of the similar 
symmetry of the RKKY coupling integral in ultrathin 
film (note that the symmetry properties of a RKKY cou- 
pling for a thin film have been mentioned first in the 
work of Wojtczak 32 ). Let us remember that the summa- 
tion over coordination zones within the atomic plane can 
be performed using the two sets of coordination radii and 
coordination numbers different for odd and even values 
of lu. However, for low values of charge carrier concen- 
tration the inverse of the Fermi wavevector is quite large 
and we are in the quasi-continuous regime so that the de- 
tailed lattice structure in the atomic plane is not of great 
importance. Thus, the symmetry £ (lu) — £ [n — lu) is 
exact for every value of p only for n even, while it is 
approximate for low p for n odd. 

In view of the experimental data, it might be instruc- 
tive to mention that the constant C/a 2 approximately 
amounts to 1.6 mJ/m 2 for (Ga,Mn)As. 

It is visible that for n = 5 to 9, the interplanar cou- 
plings are ferromagnetic (F) for the whole range of dis- 
tances (with an exclusion of wery weak antiferromagnetic 
(AF) values for n — 9, near the centre of the film). 
When n increases from 5 to 9, the coupling values are 
reduced continuously but the shape of the curve remains 
essentially unmodified. When n is changed from 9 to 
10, a QSE is visible (compare Fig. (2). Thus, the curve 
shape changes discontinuously with two minimas, sepa- 
rated by a central maximum. The coupling values are 
at the same time shifted towards stronger ferromagnetic 
values. Then, further increase of n causes the gradual 
lowering of couplings and for n = 16 the two minima be- 
come AF in character. They reach the maximum depth 
for n = 20. Switching to n — 21, the next manifesta- 
tion of QSE is visible as the curve develops an additional 
central minimum between two maximas. Once more, the 
coupling values became shifted towards F, and the two 
AF minimas become shallower. When increasing n up to 
n = 32, we notice the same tendency as before, namely 
the couplings are lowered, the central minimum becomes 
AF in character and the AF couplings near the two other 
minimas tend to rise. 

It is instructive to compare the values of interplanar 
energy obtained by performing the summation in the for- 
mula @ for the thin-film RKKY coupling integral (given 
by ©) with the some summation for the ordinary 3D 
RKKY coupling (with the same nonzero width of con- 
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tact potential)^ The results are presented in Fig. [4] for 
p = x/n, n = 20 and x — 0.05. The full circles corre- 
spond to the results derived from thin-film RKKY cou- 
pling, while the empty circles depict the interplanar cou- 
pling for 3D RKKY integral. It is visible that for the 
distances smaller that a half of film thickness, the 3D 
coupling values are shifted towards F coupling, and also 
some phase shift is present. The curve for 3D RKKY 
does not possess the symmetry with respect to the cen- 
tre of the film so that it vanishes quite fast for z/d > 0.5, 
unlike the proper thin-film RKKY coupling. In the inset 
of Fig. the values of the exchange integral itself for 
r = for ultrathin film (solid line) and bulk 3D (dashed 
line) are plotted as a function of the distance z between 
the interacting spins. Here one can easily observe the 
phase shift and the difference in magnitudes between the 
two curves. 

The interlayer coupling energy for til — 1, calculated 
on the basis of formulas (|7I9[) is a key point in our study. 
Knowing the interplanar coupling energies as a function 
of concentration of carriers and geometry of the ultrathin 
film, it is possible to search for the sets of parameters 
which result in the most robust antiferromagnetic inter- 
layer coupling. As can be seen from the formula ([9]), the 
interlayer energy is calculated by summing the interpla- 
nar coupling energies for 2ul — 1 interplanar distances, 
ranging from (77,5 + 1) d up to (2ul + ns — 1) d, with ap- 
propriate weights. Therefore, if the maximum number of 
£ (z) values, starting from z = (ns + 1) d is antiferromag- 
netic, the strongest antiferromagnetism can be expected. 
On the other hand, since the coupling for z — (ns + 1) d 
enters the sum in ([9]) with the largest weight n^ , it is also 
advisable to select ns in order to have this AF coupling 
considerably strong. For example, from the analysis of 
the Fig. [31 prepared for x — 0.05, we can expect to ob- 
tain the AF interplanar coupling for the first minimum 
when we select ns = 3 for n — 20. For thicker films, 
11 = 32, the optimum value seems to be ns = 4 (the 
first minimum) and ns = 14 (minimum in the middle, 
noticing that for ns = 13 the most important coupling 
£ (( n S + 1) d) would be only slightly antiferromagnetic). 
For the larger number of subsequent antiferromagnetic 
values for £ (z), we expect the stability of AF interlayer 
coupling for thicker magnetic layers. This thickness is 
of importance for two reasons. For one thing, with a di- 
luted system in mind, the larger n^ means larger value of 
total magnetization; for another, maximizing the critical 
temperature is also advantageous. On the other hand, it 
appears that the increase of the concentration of charge 
carriers in the system would cause the antiferromagnetic 
behaviour to diminish. This is due to the fact that the 
distance between the subsequent zeros of the function 
£ (z) decreases when p increases, so that the function 
oscillates more rapidly and the chance of having AF cou- 
plings for a few subsequent interplanar distances z de- 
creases. Let us note that the last observation seems to 
stand in apparent contradiction to the general idea that 
for small concentration of charge carriers ferromagnetism 
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FIG. 6: Interlayer coupling energy for the system character- 
ized by x — 0.05, n = 20 and ns = 3, as a function of mag- 
netic layer thickness. The two charge carriers concentrations: 
p — x/n (circles) and pLx/n (squares) are assumed. 
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FIG. 7: External magnetic field required to switch the magne- 
tizations of the layers from AF to F state as a function of layer 
thickness. The results are presented for T — » and selected 
parameters, where AF interlayer coupling is predicted. 



is favoured (because the first zero of the exchange integral 
is shifted towards larger distances and more coordination 
zones belong to this ferromagnetic range). 

This effect is illustrated in the Fig. where the in- 
terplanar couplings are plotted for three charge carrier 
concentrations: p = x/n, p = Ax/n and p = 8x/n, for 
n = 20 and x — 0.05. It might be useful to mention 
here that these values for (Ga,Mn)As would correspond 
to the (physically relevant) hole concentrations of 1.1, 
4.4 and 8.9T0 20 cm" 3 , respectively. It is clearly visible 
that for the lowest concentration of charge carriers, the 
antiferromagnetic minima are quite remarkable. When p 
is increased fourfold or eightfold, the oscillations of cou- 
pling versus interplanar distance become much faster, so 
that only two or one coupling is of distinct AF character. 
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Fig. [6] illustrates the magnetic layer thickness depen- 
dence of the interlayer coupling energy for the already 
mentioned case of x — 0.05, n = 20 and ns = 3. The 
coupling energy is plotted for two different carrier concen- 
tration regimes. The circles correspond to the constant 
carrier concentration of p = x/n, independent on the 
thickness of the magnetic layers ul- The second result, 
marked by squares, is for the carrier concentration p^x/n 
proportional to the number of magnetically doped atomic 
planes. It can be seen that for the second case, the initial 
antiferromagnetic interlayer coupling energy is lost soon 
when ri,£ increases. Although the changes are strongly 
non-monotonic, yet the AF coupling is never restored. 
The situation is very different for the constant charge 
carrier concentration. The AF coupling persists up to 
7i l = 4 and reaches the maximum strength for nj, = 3. 
Further increase in magnetic layer thickness also leads to 
change to F interlayer coupling. The similar behaviour 
can be observed for another sets of parameters regarded 
as beneficial in terms of assuring persistent AF interlayer 
coupling. 

Fig. [7J presents the critical external field B c = poH c 
required to switch the direction of layer magnetizations 
from antiferromagnetic to ferromagnetic with tempera- 
ture T — ► 0. The values have been calculated from the 
formula (fl"7|) for the parameters predicting the strongest 
antiferromagnetic interlayer coupling. In order to pro- 
vide some reference point, we calculated the value of a 
normalization constant Bq (Tl8|) for a representative DMS, 
(Ga,Mn)As, which equals Bo — 76 T, so that the critical 
fields would lie in the range of hundreds of milliteslas. It 
is visible that the critical field for each case considered is 
largest for the magnetic layer thickness or 3. If the 

lowest value of critical field is needed, together with the 
highest Neel temperature, then the best choice seems to 
be the maximum thickness of the magnetic layer which 
still provides AF interlayer coupling. 

In connection with the predictions of Fig. UJ we would 
like to mention that an attempt to swith the magnetiza- 
tion direction in antiferromagnetically-coupled multilayer 
of GaAs/(Ga,Mn)As has been made by Chung et a/.— 
In their experiment, the field necessary to reverse the 
magnetization has been estimated at about 10 mT (judg- 
ing from the magnetization curve and neutron reflectivity 
data); however, the full ferromagnetic alignment has not 
been achieved until 100 mT. It is worth noticing that in 
the experiment of Chung et al. the spacer thickness was 
12 ML, and each magnetic layer was 25 ML thick. 

In the studies of critical temperature of the system, it 
was necessary to accept a specific spin value. Therefore, 
we selected S — 5/2, which is relevant for DMS. The be- 
haviour of critical temperature of the system for the pa- 
rameters leading to AF interlayer coupling is presented 
in Fig.[5Ja). There, we plot the values of the critical tem- 
peratures for x = 0.05, n — 32 and ns = 14. Like in the 
Fig. [U the circles correspond to constant charge carriers 
density p = x/n, while the squares are for p^x/n. The 
empty symbols indicate the Neel temperatures (what cor- 



1-2 i — r 




I I I I I I I I I L 

123456789 




I i i i i i i i i i i i I 

1 23456789 10 11 

n L 

FIG. 8: Critical temperatures in MFA approximation for an 
ultrathin film with two magnetic layers. In (a) and (b) the 
filled symbols denote the Curie temperatures, while the empty 
ones are for the Neel temperatures, (a) Comparison of the re- 
sults for charge carrier concentration pLx/n (squares) and 
p — x/n (circles), (b) Results for p — x/n and various 
parameter ranges which are advantageous for AF interlayer 
coupling, (c) Comparison of the results obtained from the 
approximate Eq.(15) (empty symbols) and without this ap- 
proximation (filled symbols). 
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FIG. 9: Spatial distribution of magnetizations in subsequent 
monolayers, for various thicknesses til and for two sets of 
parameters: (a) x = 0.05, n — 32 and ns = 14, and (b) x — 
0.025, n = 28 and ns = 5. The temperature is T = (2/3) T c . 



responds to AF interlayer coupling), while the full sym- 
bols are for Curie temperature (F coupling). It is visible 
that for fixed p, the critical temperature tends to saturate 
when increasing and reaches approximately constant 
value for the magnetic layer consisting of several ML. It 
is seen that the Neel temperature for n^ — 4 is very close 
to that saturated value. For p^x/n, the critical temper- 
atures are mostly Curie temperatures and rise relatively 
fast when the magnetic layer thickness increases. 

In order to investigate more carefully the critical tem- 
peratures for the regime of fixed charge carriers concen- 
tration p = x/n, we perform the calculations for some 
other sets of parameters advantageous to AF coupling. 
The results are presented in Fig. [Sfb) . The same ten- 
dency for the critical temperature to saturate for thick 
magnetic layers is observed, and the maximum values of 
Neel temperature lie close to that limit of saturation. 

Let us mention that a similar behaviour of Curie 
temperature for thin films with RKKY interaction has 
been found in the Monte-Carlo based studies of Boselli 
et a/., 33 for fixed charge carriers concentration. How- 
ever, they considered a thin film with magnetic impu- 



rities distributed randomly in all the atomic planes and 
their study was focused mainly on non-collinear magnetic 
phases, resembling canted ferromagnetism. 

Fig. ISJc) compares of the critical temperature calcu- 
lated from the "exact" MFA formula (fT"3| and from the 
approximate formula (| 1 5|) . as dependent on magnetic 
layer thickness, for the set of parameters x — 0.025, 
n = 28 and ns — 5. It can be concluded that the as- 
sumption of uniform magnetization distribution inside 
each magnetic layer does not lead to a noticeable change 
in the critical temperature for the layers up to 4 ML 
thick. However, for thicker layers, the critical tempera- 
ture is underestimated by an approximate formula in the 
order of 5 %, which conclusion holds also for other sets 
of parameters. 

It follows from the numerical calculations that the in- 
tralayer energy is much larger than the absolute value of 
the interlayer coupling. As a consequence, it is clearly 
visible from Eq. (fT"5"| that the main contribution to the 
critical temperature originates from £i n t ra - 

As mentioned previously, in the ordered state, each 
atomic plane of the magnetic layer has its own magneti- 
zation value and some non-uniform magnetization distri- 
bution exists across the layer thickness. To illustrate this 
behaviour, we calculated such magnetization profiles for 
various thicknesses of the magnetic layers ul for the cases 
of x = 0.05, n = 32 and ns = 14 (Fig. [5)Ja)) as well as 
for x = 0.025, n = 28 and n s = 5 (Fig. Mb)). The plots 
present a magnetization distribution in a single magnetic 
layer. Each profile is calculated for the temperature of 
|T C for the given value of Ul, versus the number of the 
atomic plane counted starting from the plane situated 
closest to the undoped spacer. For the second magnetic 
layer, the values of to are either the same (for F interlayer 
coupling) or they are of opposite sign (for AF coupling) , 
so that the profile is either symmetric or antisymmetric 
with respect to the centre of the undoped spacer. It is 
noticeable that for small thicknesses, nj, = 2, the distri- 
bution is almost uniform, while the increasing thickness 
of the magnetic layer causes the profile to reshape. The 
maximum of magnetization occurs for the atomic plane 
in the middle of the magnetic layer, while the smallest 
values are reached for the planes which are the most dis- 
tant or the closest to the nonmagnetic spacer. For n/, = 8 
in Fig. M, a ) or nj, = 10 Fig.[5fb), the distribution is ap- 
proximately symmetric with respect to the middle atomic 
plane in the layer. It is noticeable that the differences 
in magnetizations for various atomic planes are strongly 
evident, since the magnetization at the boundaries of 
the magnetic layer is only a half of the maximum value 
reached in its centre. Interestingly, the profile losts its 
symmetry when we reach the maximum of possible mag- 
netic layer thickness, Ul = 9 for (a) and Ul = 11 for (b). 
Then the magnetization value in the atomic planes in the 
vicinity of the spacer is enhanced in comparison with the 
layers at the opposite side of the magnetic layer. Such an 
effect is more striking for the set of parameters x = 0.05, 
n = 32, ns = 14, where the distribution of magnetiza- 
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tion switches from almost symmetric shape with respect 
to the middle plane of the layer (for = 8) to almost 
flat one with some drop at high distances from the spacer 
(for til = 9). By approaching the critical temperature 
the profiles flatten and their magnitude is vanishing. 

IV. CONCLUSION 

The analysis of a model magnetic bilayer with ultra- 
thin film-modified RKKY interaction revealed an impor- 
tance of QSE for obtaining antiferromagnetic interlayer 
coupling. Some ranges of parameters, advantageous for 
AF were identified, with a general remark that the car- 
rier density should not be so high to make the oscilla- 
tion period excessively short. The coupling is possible 
even for relatively thick spacers (14 ML for example); 
however, the maximum thickness of each magnetic layer 
should not exceed a few ML. The critical temperatures 
were found to saturate when increasing magnetic layer 
thickness. 

The general conclusion which can be drawn from the 
study is that the occurrence of AF coupling presents a 
quite unique situation in comparison with the most com- 
mon F ordering. Thus, in order to observe the AF con- 
figuration the special set of parameters has to be cho- 
sen including the impurity concentration, its distribution 
within the thin film, the thickness and the carriers den- 
sity. Studies of sensitivity of the obtained configuration 
to the external magnetic field might be useful for practi- 
cal realization of the switching function between the AF 
and F states. On the other hand, taking into consider- 
ation the importance of selecting the spacer thickness, 
some analysis embracing the interface roughness might 
also be interesting. 

Although the magnetization profiles in thin films were 
already studied for NN interactions a long time ago 8 i 34 ' 35 , 
as far as we know for the long-range RKKY interactions 



the above presented profiles are found for the first time. 
Prediction of the profiles would need some experimen- 
tal confirmation in the future studies of DMS systems. 
Let us observe that the presence of non-uniform magne- 
tization profile causes a noticeable decrease in the total 
magnetization of each magnetic layer at finite tempera- 
tures (with respect to the prediction of the usual, bulk 
Brillouin function). 

In an ultrathin film, the characteristic features of 
RKKY exchange integral depend crucially on the 
standing-wave form of the free-carrier wavefunctions. 
They originate from the interference of waves propa- 
gating perpendicularly to the thin film plane. Thus, 
long enough phase coherence length for free charge car- 
riers is required for this picture to be valid. Taking 
(Ga,Mn)As as an example, the coherence lengths of the 
order of 100 nm have been reported in quantum wires and 
rings for millikelvin temperatures (in presence of mag- 
netic impurities)^ even though a mean free path for the 
carriers is shorter, down to an order of a lattice constant 
in the metallic-like regime^ Such a value of phase coher- 
ence length exceeds the thicknesses of the films presented 
in this paper. 

In the present work it has been assumed that the (ran- 
dom) positions of the magnetic impurity ions are uncor- 
related. However, it could be also possible to investigate 
the influence of correlated positional disorder on the mag- 
netic properties of the thin-film systems by considering 
the non-zero Warren-Cowley parameter. 
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